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MOUNTAIN BUILDING OF SOLID QUARK STARS 
Haifeng Yang 1 and Renxin Xu 1 
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ABSTRACT 

One of the key differences between normal neutron and (bare) quark stars is relevant to the fact 
that the former are gravitationally bound while the latter self-confined unless their masses approach 
the maximum mass. This difference results in the possibility that quark stars could be very low 
massive whereas neutron stars cannot. Mountains could also be build on quark stars if realistic cold 
quark matter is in a solid state, and an alternative estimation of the mountain building is present. As 
spinning compact objects with non-axisymmetric mass distribution will radiate gravitational waves, 
the equations of states of pulsars could be constraint by the amplitude of gravitational waves be- 
ing dependent on the heights of mountains. We then estimate the m aximum mou ntains and thus 
quadrupole moment on solid quark stars, to be consistent with that bv lOwenl (|2005l ) if the breaking 
strain is ~ 10 _1 , addressing that a solid quark star with mass <^ 10~ 2 M Q could be "potato-like". 
We estimate the gravitational wave amplitude of solid quark stars with realistic mountains to be order 
of h ~ 1CT 27 . 

Subject headings: dense matter — gravitational waves — pulsars: general — stars: neutron — ele- 
mentary particles 



1. INTRODUCTION 

Neutron stars, with mass of ~ M Q and radius of 
~ 10 km and thus supra- nuclear density, provide as 
unique astro-laboratory for testing fundamental theo- 
ries. Due to their high mass to radius ratio, neu- 
tron stars with strong gravity could radiate gravitational 
wav es (GWs) in a few potential wa ys: inspiralling bina- 
ries (ISathvaprakash fc Schutzll2009| ) , supernova core col- 
lapse (|Newll2003HOttll2009H . rotating deformed stars and 
GW radiation induced from o scilla tions and instabil ities 
(jAndersson fc Kokkotaslll998ft . See lAnderssonl (|2010t ) for 
a review of these topics. 

This work is focus on the gravitational waves due to the 
rotating deformed stars. It is well-known that neutron 
stars with non-axisymmetric distribution of mass may 
radiate GWs because of the time-varying quadrupole 
moment they have. The frequency of the gravitational 
waves will be twice of rotation frequency of neutron stars. 
As the the frequencies are of 10 2 ~ 3 Hz, the pulsars are 
good sources of gravitati onal waves for inte rferometers 
detectors, such as LIGO ([Abbott et al.ll2009l ) and Virgo 
(jAcernese et al.1 120081 ). Although LIGO hasn't detect 
the signal from pulsars yet, a limitation of the ampli- 
tude of the gravitational wave s has been put forward 
(|Abbott et al.ll2005l 120071 12010I ). 

One of the key motivations to detect GWs from pul- 
sar is to provide constraints on the equation of state 
of cold matter at supra-nuclear density, a challenge in 
understanding the fundamental strong interaction be- 
tween quarks. Pulsars could be neutr on or quark stars, 
even in a solid state (|Xull2003l f2009h . Can the LIGO 
detection be understood in the regime of solid quark 
stars? What about GWs if pulsars are really solid quark 
star? We are concerne d abou t these issues in this pa- 
per. lUshomirskv et all (|2000l ) first consider the elastic 
deforming neutron stars. Their estimation of the max- 
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imum quadrupole moment of the normal neutro n stars 
is 2.4 x 10 38 g-cm 2 , which has been updated bv lOwenl 
(2005) via redefining the shear modulus. They also give 
a general equation to calculate the quadrupole moment 
induced by shear modulus which ca n be u s ed to other 
stars such as solid quark stars, and lOwenl (|2005T ) esti- 
mates the quadrupole moment of the solid quark stars 
to be 2.8 x 10 41 g-cm 2 . The quadrupole moment is also 
calculated in other stellar models, such as Hybrid crys- 
talline clour-superconducting star (IKnippel fc Sedrakianl 
120091: lHaskell et aTll2007l: 1173120071) a nd Hybrid an d me- 
son condensate stars (jOwenl 120051 ) . iPitkinl (|2011l ) pro- 
vides a good review on these. 

However, all these studies assume that the equilibrium 
structure of the neutron stars is almost spherically sym- 
metric. Certainly, some solidified mountains would also 
be possible on the surfaces of sold quark stars. Unlike 
shear- modulus-induced mountains (elastic mountains), 
we suggest here latent-heat-induced mountains (solidi- 
fied mountains) there. Under this condition, instead of 
the shear modulus and the breaking strain, latent heat 
perform an important role that determine the height of 
the mountains. From this point of view, our model is 
only one parameter (i.e., latent heat) dependent, rather 
than than two parameter (i.e., shear modulus and break- 
ing strain) dependent in case of elastic mountains. 

After estimating the maximum height of the moun- 
tains, we find that if the mass of solid quark stars is 
small enough that the height is of the same order with 
the radius of the star, the star may be "potato-like" . The 
critical mass for potato- like solid quark stars is estimated 
to be ~ 10~ 3 ~~ 2 M qj which ag rees with the mutually in- 
dependent study of Xu (2010). Also it is found that, in 
reality, the actual GW amplitude of a pulsar would be 
too small to be detected with LIGO now since the max- 
imum GW amplitude requires a star to deform into a 
particular distribution of mountains. 

2. GENERAL ARGUMENTS 
Mountains on the Earth as an analogy. Before dis- 
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cussing the latent-heat-induced mountains on the solid 
quark stars, we present an estimation of the maximum 
height of mountains on the Earth using the same method. 
Consider a mountain with the height H and sectional 
area A. The earth is assumed to have mass of M and 
radius R. Local acceleration of gravity is g. If the moun- 
tain gets 5 height down, the material at the foot of the 
mountain will get energy of ASpg. While this energy is 
enough for this amount of material to melt, the moun- 
tain can no longer get higher because the mountain may 
tend to melt rather than to keep the height if it is higher 
than this. Certainly, this process of melting is somewhat 
a virtual one. We can also understand this in another 
way that this method is comparing energy between two: 
a mountain with height H and the other with height 
H — S. The maximum height of the mountain, H m , will 
then guarantee an equation of ASpgH m — ASpANa/ p, 
where AN a is the latent heat per mol (with Na the Avo- 
gadro constant) and p is the mass per mol. Therefore we 
have H m = ANA/(pg)- On the earth the material that 
construct mountains is Si02, which has p = 60 g/mol 
and AN a = 854 kJ/mol. The estimate of the maximum 
height of mountains on the earth is thus 14.5 km. The 
highest mountain known on the Earth is Qomolangma 
with the height of 8.8 km, which is the same order of the 
estimated maximum height. 

Mountains on solid quark stars. In this case, we can 
also have the equation of 
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H m A5nxniQ = ASnA, 



where n is the number density of the quark clusters, \ 
is the number of quarks in a quark cluster, mo = 300 
MeV is an assumed mass of every dresse d quark, and A 
is the melting energy per quark cluster. iMichell (1988) 
proposed that quark-alpha particle with 18 quarks is pos- 
sible to exist, while the proton and neutron are formed 
with 3 quarks. We then suggest that the x varies from 
3 to 1 8 and choose x = 10 as fiducial value. Ilshii et akl 
(2007) shows that the nucleon-nucleon potential could be 
as depth as ~ 100 MeV. The latent heat should thus be 
considered to vary from 1 MeV to 10 MeV, and A = 5 
MeV is chosen as the fiducial value for furthur calcula- 
tion. The definition of A, H m , M and R is the same as 
above. One can then has 
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It seems that there are two parameters, A and Xt but it is 
evident that the maximum height is only one parameter 
dependent, A/%, the average latent heat per quark. 

Certainly the stars are not spherical because of many 
mountains there. Stars could only be approximately 
spherical if H m <C R, otherwise, if H m /R ~ O(l), those 
gravity-free stars might be "potato-like" . For a low mass 
solid quark star with then approximately constant den- 



sity, p~4x 10 14 g/cm 3 , one has from Eq.([T]) 
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When R ~ V0.013 x 10 km = 1.15 km, we will have 
H m /R ~ O(l), and the corresponding mass for potato- 
like solid quark stars would be be M <^ (\/0.013) 3 x 
E4Mq = 0.002M Q . This result agrees with that of IXul 
( 2010), who addressed that quark stars with masses as 
low as 10~ 2 ~ 10 _3 M Q are gravity-free. An important 
astrophysical consequence of potato-like solid stars would 
be of precession with high amplitude, and we expect to 
test this by future observations. Possible low-mass com- 
pact stars is a direct consequence of the suggestion that 
pulsars could be quark stars, which could be the central 
remnants left by the deto nation of the ac cretion- induced 
collapse of white dwarfs (|Yu fc Xull201ll ). 

3. GRAVITATIONAL WAVE RADIATION FROM SOLID 
QUARK STARS 

As there are mountains on the stars, these mountains 
may attribute to the quadrupole moment of the stars and 
thus influence the amplitude of the gravitational waves 
the stars radiate. However, if a star accrete mass axisym- 
metrically, it will have an axisymmetrical distribution of 
mountains. Thus, the star will not radiate gravitational 
waves at all. It is difficult for us to know the real distribu- 
tion of mountains and amplitude of gravitational waves 
the star radiate, but we can provide an upper bound 
of the amplitude, by assuming a special distribution of 
mountains on the star. The first term that attribute to 
the quadrupole moment is Y22, so we can assume the 
stars having a height distribution of the mountains as 
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H m Re{Y 22 (0,0)} 
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where N = 15/32tt. We give the factor 1/(2VN) here to 
ensure the difference between the maximum value and 
minimum value is H ^. W e can define Q22 in the way 
I ([20M L 



of lUshomirskv et al.1 ()2000[ ). Q22 = J 5p22{r)r 4: dr. In our 
simulation, we simply define a density perturbation of 
Sp to be of delta function (only to be non zero when 
r = R). By comparing this with our distribution of the 
mountains, we can obtain Q22 = pH m R 4 /2. We have 
thus 
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Let's compare this result with that of lOwenl (|2005[ ). who 
proposed a solid quark star to have a quadrupole mo- 
ment of 2.8 x 10 41 g-cm 2 by suggesting a breaking strain 
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of 10~ 2 . Our estimation is one order higher than that 
of lOwenl (j2005) , and shows that the solidified mountains 
can provide a bigger quadrupo l e mo ment than elastic 
mountains do. iHorowitz et al.l (|2009t ) have done some 
work to show that the breaking strain on normal neu- 
tron stars can be as big as 10 -1 . Our result about the 
maximum quadrupole moment on the solid quark stars 
would agree with Owen's if the breaking strain on solid 
quark stars can also be as high as 1CP 1 . It is evident that 
both the dependence on radius and mass are the same in 
this work and that of lOwenl (|2005l) . 

What can we constrain the equation of state by 
the GW observations? To compare the theoret- 
ical quadrupole m oment with the LIGO S5 data 
(| Abbott et al.l [20Toh . we are to provide the relationship 
between Q22 and GW amplitude ho- Assuming a density 
pertu rbation of dp = Re[/r?22^22(^, </>)] ()Ushomirskv et all 
l2000f ) and following the definition of Q22 above, one 
comes to 
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where a factor l/y/2 before the amplitude of the trace- 
reversed perturbation is added for ho estimation. Similar 
treatments can also be found in oth er gravitational wave 
literatures (jJaranowski et al.lll998l ). The only difference 
is applying ellipticity e rather than quadrupole moment 
Q22- Note the relationship e = \/8n/15Q22/Izz, it is 
easy to prove that they are equivalent. 

With our estimation of the maximum quadrupole mo- 
ment into the equation above, we can obtain the maxi- 
mum amplitude of the gravitational waves, 
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LIGO is focusing on the direct detection of gravitational 
waves. Recent LIGO S5 data has g iven an upper boun d 
of the ho for 116 known pulsars ([ Abbott et al.l [20~Toh . 
The value of the upper bound of the 116 pulsars varies 
from 1CP 26 to 10~ 25 . However, from Eq.©, it seems 
that solid quark stars can have mountains high enough 
to radiate gravitational waves with amplitudes as high 
as 1CP 22 . Does this mean that these 116 pulsars cannot 
be solid quark stars? 

Should we believe that the mountains on solid quark 
stars have the maximum heights? In addition, in our esti- 
mation of Eq.©, we assume a distribution of Y22, which 
can attribute most to the gravitational waves. It is surely 
quite strange if natural mountains on stars have the ex- 
act distribution of Y22 in order to produce the maximum 
gravitational waves. 

The orogeny in the Earth's crust is powered by the 
mantle convection and thus the engagement of tectonic 
plates. What could be the force to build mountains 
on solid quark stars? Elastic energy develops when a 
star evolves, and both bulk-invariable and bulk-variable 
forces can result i n decreases of mo ment of inertia dur- 
ing a star quake (|Peng fc Xul 1200 8). This force would 



be responsible for mountain building too. We can then 
estimate the real heights of mountains from the glitch 
phenomena of pulsars below. We think glitches oc- 
cur if the momen t of inertia chang e suddenly inside a 
solid quark stars (|Zhou et alJ 12004) . The pulsars' an- 
gular momentum, under the assumption of sphere, is 
~ 3MR 2 fl/5. Making a differential of it, we have 
\6R/R\ - 5 x I0~ 7 [6(n/n)/10- 6 }. Although this con- 
sideration could not be the real heights of mountains on 
pulsars, we think that actually the height of mountains 
might be of the same order of SR. Therefore, replac- 
ing 5R with H m and inserting this into the estimation of 
Eq.((4|), we have the estimation of amplitude ho below, 
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We draw Fig. 1 to show how ho varies as the mass 
of soli d quark stars changes in the model of lLai fc Xul 
(2009), for different glitch amplitude SQ/Cl. In the calcu- 
lations, we choose the numbers of quarks in a quark clus- 
ters to be 18 and the potential Uq to be 50 MeV. We can 
see that if our estimation of the real height of mountains 
is valid, it would be natural that LIGO still hasn't detect 
the gravitational waves directly, since the maximum GW 
amplitude presented in Eq.([5]) requires maximum height 
and a particular Y22-distribution of mountains. What if 
the mountains haven't the Y"22-distribution? This will be 
discussed in the next section. 

4. GRAVITATIONAL WAVES INDUCED FROM RANDOM 
DISTRIBUTION 

Above we give an approach to the actual amplitude of 
GWs using the glitch phenomenon. However, the ampli- 
tude of GWs from quark stars depends also on another 
factor. Besides the maximum height of mountains, the 
distribution of mountains plays a key role as well. In the 
previous section, we give an estimation of the maximum 
amplitude of GWs by assuming a specific distribution 
described in Eq. ((2]). Such a hypothesis is too strong to 
approach the physical circumstances. In this section we 
will discard this assumption and consider another likely 
distribution. 

It is certainly very difficult to know the real distribu- 
tion of mountains on a solid quark star. Nevertheless, an 
idea comes out that we can give a random distribution 
to approach the actual situation. Yet another problem of 
the definition of random distribution emerges. The most 
strict and physical definition should satisfy the following 
requirements. (1) The height of mountains should varies 
from to the maximum height of mountains. (2) The 
surface of the star should be continuous, say the func- 
tion H{9, (f>) which describe the distribution of mountains 
should be infinitely differentiable. (3) The third require- 
ment comes not from mathematics, but from physics: the 
partial derivative of H(6, (j>) should not be too high, oth- 
erwise a mountain which is quite precipitous would tend 
to fall down. 

With all the considerations above, we find it is impos- 
sible or at least very hard to use the most strict definition 
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of random distribution to approach the real amplitude of 
GWs. Nonetheless, we would have a order-of-magnitude 
estimation of the distribution. The only thing we do care 
is the quadrupole moment the distribution induces. We 
then use another definition which is not that strict but 
we expect the latter definition can give a similar predic- 
tion or at least in the same order. We simply divide our 
star into some small pieces, say 100 x 100 for example. 
Each piece has height varying from zero to the maximum 
height of mountains in a uniformly random way. This is 
the definition of H r (8, <fi) which is discussed below. 

As a scatter function defined on sphere, the function of 
H r (9,(f)) can be expanded using the spherical harmonics 
as following, 

I . I a 

Each mode contributes to corresponding mode of GWs. 
The mode with the sphere distribution of 122(0, 0) is be- 
lieved to contribute GWs with the biggest amplitude. 
To pick up this mode we use the orthogonal property 
of the sphere harmonics. We multiply the H r (8,4>) with 
Re[F22(0, </>)] and calculate the integral on the sphere. 
The result is regarded as the factor before Y22 in Eq. 
([2]), and we can calculate the quadrupole moment ac- 
cordingly. 

We divide the star into 100 x 100 pieces, calculate 10000 
times in the simulation, take the average of the results 
and finally obtain the estimation of the quadrupole mo- 
ment of the order 10 38 g • cm 2 . Correspondingly, the am- 
plitude of GWs is of the order 10~ 27 . Considering the 
upper bond given by LIGO S5 data varies from 10 -26 to 
10~ 25 , our results would be good enough to explain the 
conflict between maximum estimation and the observa- 
tions. 

5. CONCLUSIONS AND DISCUSSIONS 

We present a new method of building mountains on 
solid quark stars. It is induced by the latent heat rather 
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Fig. 1. — Estimation of amplitude of gravitational waves from solid quark stars in lLai fc Xul (200$ ) model (N q = 18 and Uq = 50 MeV). 
This figure shows how GW amplitude, ho, varies while the mass of stars changes. The line on the top is ho estimated from the maximum 
height of mountains with the distribution of Y?,2- The five lines under it represents for the hg estimated from the consideration of glitch. 
See the text for details. 



